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ABSTRACT 

Thanks to their past history on the main sequence phase, supergiant massive stars 
develop a convective shell around the helium core. This intermediate convective zone 
(ICZ) plays an essential role in governing which g-modes are excited. Indeed a strong 
radiative damping occurs in the high density radiative core but the ICZ acts as a 
barrier preventing the propagation of some g-modes into the core. These g-modes can 
thus be excited in supergiant stars by the K-mechanism in the superficial layers due 
to the opacity bump of iron, at logT = 5.2. However massive stars are submitted to 
various complex phenomena such as rotation, magnetic fields, semiconvection, mass 
loss, overshooting. Each of these phenomena exerts a significant effect on the evolu- 
tion and some of them could prevent the onset of the convective zone. We develop a 
numerical method which allows us to select the reflected, thus the potentially excited, 
modes only. We study different cases in order to show that mass loss and overshooting, 
in a large enough amount, reduce the extent of the ICZ and are unfavourable to the 
excitation of g-modes. 
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1 INTRODUCTION 

Hipparcos data set has provided new samples of period- 
ically variable B-type stars. Among them, about 30 pe- 
riodically variab l e supergiant stars have been detected 



riodically variab l e supergiant stars have been detected 
JWaelkens et all Il998l: lAerts et all Il999l : iMathias etall 



I2OOII ). iLefever et al.l (|2006l ) have recently reexamined this 
sample and they derived the atmospheric parameters 
through line profile fitting. From the location in the (log (7, 
logTeff) diagram (Fig. [!}, they suggest that the variability 
is due to n o n-rad ial pulsations excited by the Ac-mechanism. 
ISaio et al] (|2006l ) also reported the discovery of p and g- 
mode pulsations in a B supergiant star HP 163899 (B 2 Ib/II 
iKlare fc Neckel Il977l : ISchmidt fc CarruthersI 1 19961 ) which 
has been observed by the MOST satellite. For this star 48 
frequencies have been detected (< 2.8 c/d) with amplitudes 
of a few milli-magnitudes. 

At first sight, the presence of g-mode pulsations in a su- 
pergiant star is quite unexpected. Indeed such stars present 
a rather condensed radiative helium core with a very large 
Brunt- Vaisala frequency which produces an important ra- 
diative damping. This radiative damping is so strong that 
no g-modes ente ring the core should be excited. However 
ISaio et al] (|2006l ) have shown that the presence of excited 
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g-modes is indeed possible thanks to an intermediate con- 
vective zone (ICZ) which prevents some of the modes from 
entering the radiative damping core. In that case, the k- 
mechanism in the superficial layers can be sufficient to ex- 
cite the modes. The presence of this ICZ is therefore crucial 
to have excited g-modes in supergiant stars. 

However massive stars are submitted to various com- 
plex phenomena such as rotation, strong mass loss, magnetic 
fields, semi-convection and overshooting. Each of these phe- 
nomena exerts a significant effect on the evolution and some 
could prevent the formation of an ICZ. An enlightening and 
still wid ely used review on mas sive star evolution is to be 
found in Chiosi & Maeder (1986') while a more recent review 
has been made by Maeder. (1998 ). We shall here limit our 
discussion to slow rotators, slow enough to be able to ignore 
the additional mixing induced by rotation. Magnetic fields 
will also be ignored. 

We shall here consider models on the main sequence 
(MS) or close to it and analyze the effect of mass loss and/or 
overshooting. Mass loss in duces a faster receding convective 
core during the MS phase (|Chiosi fc Maederll 19861 ) . If a large 
enough mass loss rate is taken into account during the main 
sequence, no ICZ can be formed during the post-MS phase. 
Overshooting will be discussed according to two different 
aspects: (1) overshooting during MS can prevent the for- 
mation of an ICZ during the supergiant phase; (2) with a 
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Figure 1. log g — logTeff diagram in w hich the sampl e of pe riod- 
ically variable B type supergiants from pLefever et ahl 1I2OO6I ) and 
some evolutionary tracks from 10 to SOMq computed with CLES 
are displayed. The MOST star is also sho wn with the error box de- 
rived from photometric studies ( Schmid t &: Carruthersl ll996j) by 
ISaio et al] l l2006f) . iLefever et al.l I I2006.') have divided their sample 
into 4 groups depending on the reliability of the derived stellar 
parameters (mainly T^ff and logg). Red points are the results 
considered as very reliable, restricted reliability is shown in blue. 
Green points stand for the stars for which they had no means to 
derived accurate stellar parameters, the results are thus consid- 
ered as unreliable and must be taken with caution. They used 12 
comparison stars in which 9 of them appear to be periodically 
variable too (magenta points). 



large amount of overshooting a star located in the super- 
giant region could still be on the MS, i.e. with a convective 
hydrogen-burning core. Although the location of these stars 
in the (logg, logTeff) diagram could also be attributed to 
stars undergoing central helium-burning at the blue part of 
the He loop, we shall concentrate in this paper on the H 
burning phases: either the main sequence or the supergiant 
phase. 



2 MODELS 

The models were compute d with the Code 
d'Evolution SteUaire, CLE S (IScuflaire et all 
used the OPAL opacities (llglesias fc Rogers! Il99t 



Jiegeois 
We 



pleted with the [Alexander fc FergusonI ( 19941 ) opacities at 
low temperature. We adopted the CEFF equation of state 
l|Christensen-Dalsgaard fc Daeppenlll992l ) and the old stan- 
dard heavy-element mixture (Grevesse fc Noels 1993). 
Boundary conditions at the surface are derived from a ra- 
diative gray model atmosphere computed in the Edding- 
ton approximation and with the OPAL equation of state. 
The a diabatic oscillation fre quencies are computed with 
LOSC (|Scuflaire et al.ll2007al ) and the excitation of oscilla- 
tion modes is comp uted using the non-adiabatic code MAD 
jPupret et al.ll2003l ). 



3 ROLE OF AN ICZ IN THE EXCITATION OF 
G-MODES IN SUPERGIANT STARS 

3.1 Internal structure 

B supergiant stars are post-MS massive stars which are 
burning hydrogen in a shell surrounding the helium core. 
Some of them are already burning helium in the core. The 
He core is very dense and is contracting whereas the low 
density radiative envelope is expanding. The structure of the 
star in the vicinity of the hydrogen burning shell depends 
on the structure of the star during the MS. In particular, 
depending on this past history, an intermediate convective 
zone can develop above the core. 

The formation of the ICZ requires a region where the 
neutrality of temperature gradients is reached: Vrad = Vad 
during the MS phase. This region is formed by the decrease 
of the radiative and the adiabatic temperature gradient dur- 
ing the MS. The transformation of H into He decreases the 
opacity and thus the radiative temperature gradient which, 
in less massive stars, leads to a receding convective core. 
However the larger the initial mass, the higher the central 
temperature and massive stars have therefore a large radia- 
tion pressure which increases during MS. This results in a de- 
crease of the adiabatic temperature gradient and the convec- 
tive core either recedes much more slowly or even increases in 
mass during the MS. When this happens, a semi-convective 
zone appears, and adopting the Schwarzschild criterion for 
convective stability, a partial mixing occurs in a region sur- 
rounding the convective core with a radiative temperature 
gradient very close to the adiabatic one (Fig. (21). In the post- 
MS phase, when H-shell burning starts in this region, the 
radiative gradient becomes very rapidly larger than the adi- 
abatic one since L/m becomes very large (Vrad oc nL/m) 
and an intermediate convective zone appears in, or near, 
the H-burning shell. 

However that MS 'neutral' region as well as the post- 
MS ICZ can be affected by various physical processes occur- 
ring during MS. If a large enough mass loss rate is taken 
into account during the MS phase, the central tempera- 
ture increases less quickly. Hence, the adiabatic gradient 
no longer decreases significantly. Mass loss also results in 
a convective core receding more quickly a nd does not leave 
behind it a region where Vrad — Vad (|Chiosi fc Maederl 
119861). Therefore, with significant mass loss, no ICZ appears 
during the post-MS phase. Overshooting during MS can 
also prevent the formation of an ICZ during the supergiant 
phase. Finally, the criterion used to define the boundaries of 
the convective zones: Schwarzschild Vrad = Vad or Ledoux 
Vrad = Vad + (/3/(4 — 3 /?))Vn = Vl also affe cts the appear- 
ance or not of an ICZ (|Lebreton et ahlboOSl ). 



3.2 Excitation of g-modes 

The structure of a supergiant star is similar to the one 
of a red giant: it is characterized by a high density con- 
trast between the small size core and the expanding enve- 
lope. This problem has b e en discussed for gi a nt sta rs by 
iDziembowskH (|l97ll . [l973); iDziembowski et al.l (|200ll ). The 
regions of mode trapping are shown in the propagation di- 
agram (Fig. [3} in which the dimensionless squared Brunt- 
Vaisala (A'^^) and Lamb (L|) frequencies are plotted as a 
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Figure 2. Radiative and adiabatic temperature gradients 
(dashed and solid line respectively) in a model of I6M0 on the 
MS phase. The dotted line stands for the hydrogen abundance. 
During the MS both the radiative and the adiabatic gradients de- 
crease and the region with the neutrality of gradients is formed. 



function of logT. We first consider a model without ICZ 
(dotted line). The behaviour of Le = + 1) c/r is quali- 
tatively similar from star to star, being infinite at the center 
and decreasing monotically towards the surface. This is not 
the case for the Brunt- Vaisala frequency: 



d In p 1 dlnP 



dr 



Fi dr 
V + Vu 



(1) 



for a fully ionized gas (2) 



which is affected by the presence of convective region where 
A'^^ is close to zero and by the mean molecular weight gra- 
dient, which leads to bumps in the Brunt- Vaisala fre- 
quency. Moreover, the Brunt- Vaisala frequency depends on 
g^p/P, which in supergiant stars, due to the contrast be- 
tween the contracting helium core and the expanding enve- 
lope, takes huge values within the radiative core. Hence the 
eigenmodes of any moderate frequency (u < 100) show a g- 
mode behaviour in the core, where a <^ N, Lf. all non-radial 
modes are mixed modes. Moreover, the eigenfunctions have 
a large wavenumber, k, in the core since: 



if (T < N,Lf 



(3) 



In this radiative g-mode cavity, where oscillations present 
short wavelengths, a strong radiative damping always oc- 
curs. Indeed, locally, the mechanical energy lost per unit 
length by the mode during each pulsation cycle due to ra- 
diative heat exchanges, is approximately given by: 



5Td'^{5T/T) 



crdlnT/dr T 



£{£+!) — — 



dr^ 



167rac 



(4) 



This expression will be detailed more rigorously in Sect. 14.11 
It is obtained by substituting one of the dominating term in 
Eq. |6] (last term) into Eq. [7] This radiative damping of non- 
radial modes is thus very large due to the huge values of A''^ 
in the core, and it seems that there are no hopes to observe 
such modes in post-MS stars. 

However, in the reasoning above we have assumed that 
the eigenfunction square modulus \5T /T^ is not negligible, 
which is not always the case due to the presence of an ICZ 
(Fig. m solid line). In a convective zone A*"^ < and g-modes 
are evanescent. Therefore, the ICZ can act as a potential bar- 
rier: some modes can cross it, others are reflected. This is 
illustrated in Fig. |4] For the mode that crosses the convec- 
tive barrier (gray line), the amplitudes of short wavelength 
oscillations are significant in the radiative core and strong 
radiative damping ensues (see the gray line on the work in- 
tegral, shown on Fig. [5]): this mode cannot be observed. The 
other mode, whose frequency is very close, is refiected on the 
convective barrier, i.e. has small amplitudes in the radiative 
core (black line) . In this case the radiative damping remains 
small compared to the K-mechanism occurring in the iron 
opacity bump near the surface (see the increasing work at 
logT ~ 5.2, black line on Fig. [5]); this mode is unstable 
and could be observed. In the dense spectrum of non-radial 
modes, some are reflected on the convective shell and can be 
excited and observed. In order to select the refiected modes, 
i.e. those with small amplitudes in the radiative core, we im- 
pose a rigid boundary condition at the base of the ICZ (see 
Sect. |4|. If the wavelength of the considered mode is a lot 
smaller than the width of the ICZ, this condition allows to 
achieve small amplitudes in the radiative core. It is useful to 
note that, if only reflected modes are considered, the mode 
propagation cavity exactly looks like a MS B star: because 
of the reflection, the radiative core can be somehow forgot- 
ten and the bottom of the cavity is a convective region. As a 
consequence, the frequencies of reflected modes behave like 
frequencies in a MS star. Limited to the reflected modes, the 
frequency pattern is sparse enough and asteroseismology of 
blue supergiants becomes possible. 



SUPERGIANT STARS: 
METHOD 



NUMERICAL 



The main numerical problem when solving the non-adiabatic 
pulsation equations in supergiant stars is the large number 
of mixed modes with a huge number of nodes. Indeed, for a 
given period range, the number of nodes, n, depends on the 
Brunt- Vaisala frequency: 



Pn.l 



27r2(n + I +e) 



VWTT)r^^dr 



(5) 



For a fi Cephei type mode with a period of the order of the 
fundamental radial mode, the number of nodes is of the order 
of 100, or even 1000 for SPB type g-modes. The problem is 
numerical: it is difficult to solve the differential equations 
in such a cavity with so large a number of nodes. We shall 
however see that it is possible to avoid the computation of 
the full dense spectrum of non-radial modes by preselecting 
the modes that are reflected at the bottom of the ICZ (and 
potentially excited). 

First of all, in order to select the reflected modes only. 
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Figure 3. Propagation diagram for two supergiants model 
(13M0, M = 0, logTeff = 4.26, log L/Lq = 4.50, solid line) 
and {16Mq, M = IQ-^Me/yr, logTeff = 4.26, log L/Lq = 
4.47, dotted line). The squared dimensionless Lamb frequencies 
(L| B? / (GM); similar for both models) are shown by solid lines 
for 1 = 1. The squared dimensionless Brunt- Vaisala frequencies 
(iV^ B? /{GM)) are shown by solid and dotted lines for the model 
with ICZ (without mass loss) and without ICZ (with mass loss) 
respectively. The horizontal dashed line corresponds to a typical 
SPB type g-mode. 



we impose a rigid boundary condition at the bottom of the 
convective shell: = 0. This condition leads to small am- 
plitudes of all eigenfunctions in the radiative core since \P'\ 
is much smaller than \5P\. As a consequence, \5P\ ~ pgl^rl 
stays also small. This condition allows us to solve the 'en- 
velope problem' individually, without taking care of what is 
happening below. All the computations are performed with 
the standard fully non-adiabatic code. Artificially imposing 
a boundary condition inside the star will of course not pro- 
vide the real frequencies of the full star. But as the frequency 
spectrum of the full star is very dense, there always exists 
a mode with a frequency close to the one corresponding to 
such a rigid reflection and the error introduced by this ap- 
proach is very small. 

Once the non-adiabatic pulsation equations are solved 
in the envelope, we cannot completely disregard the radia- 
tive core because reflection is only partial. In the radiative 
core, we have seen that the eigenfunctions have a huge num- 
ber of nodes (up to n « 1000!). This makes very difficult to 
accurately solve the differential problem in this region. But 
at the same time, two approximations are valid with high 
accuracy in this region. 



4.1 Quasi-adiabatic treatment 

Firstly, the internal energy of these deep layers is very 
high, which justifies the so-called quasi-adiabatic approxi- 
matio n. This approximatio n can be presented in different 
ways (jPziemb owski 1977: Unno et al.l 1989) , but leading al- 
ways to the same final result. We can present it also as 
follows. First, the adiabatic problem is solved to get the 
eigenfunctions 5P, ■ ■ ■ Next, these adiabatic eigenfunc- 
tions are used to determine 5L from the perturbed diffusion 
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Figure 4. Radial displacement eigenfunction versus the temper- 
ature of the star. The ICZ is shown by vertical solid lines. Two 
modes are shown, the gray one crosses the ICZ and enters the 
radiative core in which it has large amplitude, whereas the other 
(black dashed line) is reflected on the ICZ. In that case, the am- 
plitude in the core is small. The amplitudes in the envelope are 
roughly the same. The behaviour of the eigenfunctions in the cen- 
ter is an artefact coming from our use of the asymptotic approxi- 
mation (Eg. 1101 where and r come to zero). Top Representation 
of the whole star. Bottom Zoom on the central regions. 



equation: 



+ 3 



ST 5k 



e{£ + i) 



a , dST/dr 
dT/dr 



+ 



(6) 



Finally, with all these ingredients the work integral can be 
computed and divided by the inertia to get the damping rate 
of the modes, which gives in a pure radiative zone without 
nuclear reactions: 



V = 



Jo T \ 



asL 

dm 



e(e+i)L T' 

iTvpr'-' dT/dlni- 



dm 



2ct 



(7) 



where a (resp. 77) are the real (resp. imaginary) parts of 
the angular frequency (time-dependence: exp{iat — Tjt)). The 
quasi-adiabatic treatment is used for the integration in the 
radiative core (r < ro) and the full non-adiabatic eigenfunc- 
tions are used for the envelope (r ro). 
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the dominating terms in the asymptotic limit gives thus: 
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Figure 5. Work integral versus the temperature of the star. The 
work is positive at the surface if the mode is excited. The reflected 
mode (in black) is thus excited whereas the non-reflected one (in 
gray) suffers the strong damping in the radiative core. 



4.2 Asymptotic treatment 

Secondly, another important simplification is to use the 
asymptotic theory in the radiative core. This theory applies 
perfectly there since the wavelength of the eigenfunctions 
is by far smaller than the scale heights of different equilib- 
rium quantities. A full non-ad i abatic asymptotic treatment 
was derived by iDziembowskil (|l977l ). Here we use instead 
the standard adiabatic asymptotic theory, which gives the 
following expressions for the radial (^r) and transversal (^h) 
co mponents of the dis placement far from the edges of a cav- 



co mponents ot tne dis 
ity (|Unno et al.ll 19891 ) 




k 

r 

r 



In a g-mode cavity (cr << L 



k. 

r 

k 

r 



(8) 



(9) 



(10) 



(11) 



where the local radial wavenumber is given by kr = 
^l{£ + 1) N I (cr r). The K constant is obtained by applying 
the continuity of &P at the bottom of the ICZ (ro), which is 
equivalent to the continuity of P' because of our rigid bound- 
ary condition. We note from these equations that IC/il/lCi-l ~ 
Njo » 1 and \P'\/\dP/dr^r\ « Nar/g ^ a/Le « 1. 
Hence SP/P ~ (dlnP/dr) (the same for T and p) and 
\5p/p\ « £{£ + 1) (near incompressibility) , so that 

d^r/dr -£{£+!) ^h/r. 

Substituting Egs. llOl and llll in Eq. [6]and keeping only 



^ d (ST/T) 
T ~ dlnT 



£(£+1) k 
r 

£{£+!) f^-l) |. 



(12) 



We obtain then from the equation of energy conservation 
and neglecting the transversal component of the divergence 
of the flux (since dluT/dr and k,. > 1(1 + l)/r): 



ioTSs 



dSL 
dm 

^W^+1)]^^ /V^d 
47r V V 



p3(j5j,ll 



£{£ + !) 
4-K 



sm 

Vad 
V 



kr dr 



u2 y.3p 



(13) 



The contribution of the radiative core to the numerator 
of Eq. ([7| is thus simply given by: 



"0 ST dSL 



dm 



T dm 

'£ + !) r pg Vad - V 

Jo 



V,aN^L (k] dr 



P V 

[£(^ + l)]«/2 ro Vad-V V^iNgL 
Jo 



2cr3 

and for the denominator: 



V 



p j,5 



2a' 



CI dm 



47r K' 



Jo 



kr dr . 



dr (14) 



(15) 



Thes e equations are p erfectly co mpatible with those 
given in IVan Hoolst et aU ([1998) and iDziembowski et al.l 
(2001). Finally, it is important to emphasize that using La- 
grangian or Eulerian formalisms can lead to different nu- 
merical results. In our first computations, we used a La- 
grangian formalism, but it did not lead to the appropriate 
evanescent behaviour of the eigenfunctions in the ICZ; in- 
stead they showed a large wavelength oscillation. This comes 
from the fact that no control of the Brunt- Vaisala frequency 
is possible in a Lagrangian formalism. Because of numerical 
truncation errors, the code does not know if A'^^ is slightly 
positive or negative in the ICZ. With an Eulerian formalism 
instead, A''^ appears explicitly in the movement equation, 
and the correct value can be attributed to it. As was already 
pointed out by Dziembowski years ago (private communica- 
tion), we emphasize thus that it is much better to use an 
Eulerian formalism for the finite difference scheme inside a 
deep convective zone. 

Our numerical method give us therefore the frequen- 
cies of the refiected and potentially excited modes by solv- 
ing numerically, with the fully adiabatic code, the 'envelope 
problem' first. We then compute the damping rate of the 
modes taking into account the envelope and the core for 
which we made some valid approximations: the use of the 
quasi-adiabatic approximation and the eigenfunctions of the 
asymptotic theory. 
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5 RESULTS AND DISCUSSION 

New samples of B supergiant stars have been observed 
and it is now co mmonly admitted that they present non- 
radial pu lsations ( Waelkens et al. 'logs'; Aerts et aJj 
iMathias et al. 2001; Saio ot al. 2006; Lefevcr et alJli 
The Canadian MOST satellite observed the B supergiant 
HD 163899 in June 2005 during 37 days and detected 48 
frequencies < 2.8 c/d with amplitudes of a few milli- 
magnitudes. Blue supergiants are however characterized by 
a radiative core in which the Brunt- Vaisala frequency takes 
huge values. A strong radi a.tive damping ens ues and no such 
modes should be observed. ISaio et ahl l|2006t ) have suggested 
that the presence of an ICZ around the core prevents some 
modes from entering the core. In that case the radiative 
damping is largely reduced and the K-mechanism in the su- 
perficial layers can excite the mode. Some physical processes 
could however prevent that ICZ from developing; this is the 
case with mass loss and overshooting. We shall successively 
envision six different physical situations. 



5.1 A supergiant star without mass loss and 
without overshooting 

As mentioned, the excitation of the g-modes in supergiant 
stars is due to the opacity bump of iron in the surface 
layers thanks to the blocking effect of an ICZ (Sect. [S}. 
ISaio et akl (|2006l ) computed evolutionary models for a mass 
range of 7 ^ M/Mq ^ 20 with an initial composition 
of (X,Z) = (0. 7,0.02) and the OPAL opacities. Their models 
fit quite well the observations, though a frequency gap be- 
tween p and g-modes still appears in the theoretical excited 
frequencies compared to the observed ones. They derived 
the e stimated effective temperatu re from photometric anal- 
ysis (ISchmidt fc Carrutherill996l '): for a B2Ib/II type star, 
logTeff is roughly between 4.22 and 4.32. The model which 
best reproduces the observed frequency range is a model of 
15M0 at logTeff ~ 4.36 for which they assume that the fre- 
quency gap would be filled at least partially by rotationally 
splitted modes. We have computed models with mass from 
10 to 18M0. Fig. [S] shows evolutionary tracks with instabil- 
ity region for p and g-modes from 1 to 14 M(;) (top panel). 
In agreement with the work of Sai o et al.l we find excited 
g-modes on the supergiant phase. 



5.2 A supergiant star with mass loss and without 
overshooting 

The occurrence of stellar winds is ubiquitous for stars of al- 
most all mass and constitutes a key aspect of stellar physics. 
In terms of the mass loss integrated over the whole star life 
time, it is most extreme in massive luminous stars, where the 
winds are almost certainly radiatively driven. Most of hot 
massive stars (O, B, WR) indeed present a wind essentially 
due to the radiation pressure: winds ar e driven by the reso- 
nant scattering of the stellar radiation (jLamers fc Cassinellil 
1 19991 ). When taking mass loss into account, the central tem- 
perature and the radiation pressure increase less rapidly in 
the course of central H-burning. As a result, the adiabatic 
temperature gradient decreases less and less as the mass loss 
rate increases. If a large enough mass loss rate is assumed 
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Figure 6. Evolutionary tracks and stability data in the HR di- 
agram. Th e black box is the error box of the MOST star from 
ISaio et alJ (2006). Top Without mass loss (from 10 to I4M0). 
All sequences present excited g-modes (crosses) on the supergiant 
phase. Bottom With M = lO'^MQ/yr (from 16 to 20Mq): 
excited g-modes are present only for a mass higher or equal to 
17M0. 



during the MS, no ICZ is formed during the supergiant phase 
(Chiosi & Maeder 1986). 

We adopted a mass loss rate of lO~^M0/yr which is 
slightly larger than most observed MS mass loss rates in or- 
der to emphasize the effect of mass loss o n the excita t ion of 
g-mo des. Indeed, the mass loss recipe of lVink et al ] (|200G| . 
I2OOII ) gives mass loss rates of the order of 10 ^M0/yr to 
lO"*M0/yr during the MS and of about lO"'^M0/yr on 
the post-MS phase for 14M0 and I6M0 stars. But the real 
mass los s rates are still un der discussion until now (see for 
example Bouret et al.l 12005 



jEl : 

I2OO6I : lOskinova etalTbOOTl : 

ences therein). 

We have computed evolutionary tracks with and with- 
out mass loss. Their location in the HR diagram is shown 
on Fig. |6l tracks on top panel (resp. bottom panel) are 
computed without (resp. with) mass loss. We performed 
non-adiabatic computations to determine whether unstable 
modes were present or not. For the sequences computed 
without mass loss, there are indeed excited g-modes dur- 
ing the supergiant phase for models within the error box for 
the MOST star HD 163899. Even at IOM0 the supergiant 
phase is characterized by excited g-modes. This is true for 
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Figure 7. Frequency distribution of excited p and g-modes for su- 
pergiant models of ISMq computed without mass loss (left panel) 
and of 18M0 computed with M = IO-'^Mq/jt (right panel). The 
horizontal dotted lines stand for the observed frequency of MOST 
jSaio et al.ll2006h . 



all higher masses. On the other hand, in the sequences com- 
puted with mass loss, we find excited g-modes only for stars 
initially more massive than about 17M0, i.e. with an instan- 
taneous mass of 13.7M0 on the supergiant phase. With an 
even higher mass loss rate, this value becomes larger and 
larger. These sequences without any excited g-modes dif- 
fer by the absence of ICZ on the supergiant phase, due to 
the mass loss effect during the MS phase. Hence, all the 
modes enter the radiative core and suffer the strong radia- 
tive damping. The frequency distribution of the theoreti- 
cal excited_modes is shown in the same kind of figure as 
m laaio et al.l ()2006l ) but here comparing without and with 
mass loss (Fig. [Tjl ■ In both cases the frequency distribution 
of the theoretical excited modes is in good agreement with 
the observed frequencies. 

5.3 A supergiant star without mass loss and witli 
overshooting 

The convective stability criterion fixes the boundary of the 
convective core as the layer where the buoyancy acceleration 
vanishes. The velocity of the globule is however different 
from zero and the globule can penetrate the radiative re- 
gion. Such an enlargement of the mixed core makes the star 
more luminous, it also increases the core H burning lifetime 
while, in the HR diagram, the MS track is more extended. 
We model here overshooting by just extending the mixed 
region but keeping V = Vrad in the overshooting region. 
When assuming a large enough overshooting during the MS, 
the MS 'neutral' region (where Vrad — Vad) does not exist 
above the convective core and this can prevent the forma- 
tion of an ICZ during the post-MS phase. No g-modes are 
thus found to be excited. We have checked the presence of an 
ICZ by progressively increasing the overshooting parameter 
from one evolutionary sequence to the next (Fig. [8]). For a 
model of I2M0 the ICZ is well develop for a small amount 
of overshooting Oov = 0.2, it is much smaller for 0.3 and it 
completely disappears for 0.4. In the cases considered here. 
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Figure 8. Evolution of the mass extension of the ICZ during the 
supergiant phase. On the x-axis we have the effective temperature 
decreasing after the MS. Each point of this axis is a model. For 
each model the mass fraction of the limits of the convective zone 
are indicated on the y-axis. For a small amount of overshooting 
c«ov = 0.2, the ICZ is well develops, it is much smaller for 0.3 and 
it disappears for 0.4. 



the overshooting is however small enough to keep the error 
box well within the post-MS phase. 

5.4 A main sequence star without mass loss and 
with overshooting 

Another possibility to solve the problem of the presence of 
excited g-modes would be to bring MS evolutionary tracks 
into the error box of HD 163899. This could be achieved by 
including larger overshooting in MS models. Indeed, with 
overshooting, the MS phase of evolution reaches lower effec- 
tive temperatures. We computed evolutionary tracks with 
different overshooting parameter, ranging from Oov = 0.2 to 
0.5 (Fig. [9]). MS evolutionary tracks cross the error box for 
an overshooting parameter equal to or larger than 0.3. 

On the MS phase, massive stars present a convective 
core surrounded by a radiative envelope. The Brunt- Vaisala 
frequency is therefore zero in the core and no damping 
can occur. The K-mechanism due to the Fe-opacity bump 
at logT = 5.2 excites p and g-modes. We performed non- 
adiabatic computations which revealed excited g-modes in 
all the sequences. The spectrum of the theoretical excited 
modes is shown on Fig. [TD]for = 0.4 and 0.5 during the 
MS (decreasing effective temperatures) and near the turn off 
(increasing effective temperatures). The agreement in the 
mode spectrum is however not as good as it was for the 
supergiant model. 

5.5 A supergiant star in the helium burning phase 

Blue supergiants could also be He-burning stars. Their struc- 
ture shows a small convective core surrounded by a radiative 
envelope in which an ICZ may still exist. As long as such an 
ICZ is present g-modes can still be excited. After the dis- 
appearance of the ICZ, the radiative damping is no longer 
reduced by this reflective barrier and g-modes are all stable. 
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Figure 9. Main sequence evolutionary tracks computed with dif- 
ferent overshooting parameters ranging from Oov = 0.2 to 0.5. 
The black box is the error box of HD 163899. Main sequence 
evolutionary tracks with at least Qov =0.3 cross this error box. 



In a forthcoming paper we shall analyze the excitation of 
g-modes in He burning models. 

5.6 Ledoux's criterion 

The hydrogen rich layers surrounding the core during the MS 
phase are potentially unstable towards convection and a par- 
tial or semi convective mixing is supposed to take place in or- 
der to achieve neutrality versus the convective criterion. This 
criterion is either the Schwarzschild criterion or the Ledoux 
criterion and we choose here the Schwarzschild one for which 
a region is convective if Vrad Vad- However, adopting the 
Ledoux's criterion (Vrad = Vl = Vad + (/3/(4 - 3f3))\/^, 
where V^ = dln/i/dlnP and f3 is the gas pressure to total 
pressure ratio) alters the structure of the star since it adds 
to the adiabatic temperature gradient the mean molecular 
weight gradient which is formed on top of the convective 
core during MS. The MS convective core is therefore not 
altered by the choice of criterion (V^ = 0), but adopting 
Ledoux's criterion a ffects the size and th e location of the 
ICZ on the post-MS (|Lebreton et al.1 12008''). During the MS, 
the presence of the mean molecular weight gradient prevents 
the formation of a region in which the neutrality of gradi- 
ents is reached (Vrad — Vad)- However, outside the core, the 
radiative gradient increases due to a larger opacity resulting 
from the increase in the hydrogen abundance. As a result, 
an ICZ appears above the region of variable mean molecular 
weight (Fig. 1111 top panel) . This ICZ develops in an homo- 
geneous region where no nuclear reactions take place. The 
bottom of the ICZ corresponds to the base of the constant 
hydrogen rich envelope, which corresponds to the maximum 
extent of the MS convective core. On the post-MS phase this 
ICZ remains at the same location, well outside the nuclear 
burning shell fFig. 1111 bottom panel). 



6 CONCLUSIONS 

The formation of an ICZ during the post-MS phase is a re- 
sult of the MS evolution. The presence of excited g-modes in 
B supergiant stars depends therefore on physical processes 
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Figure 10. Frequency distribution of the excited p and g-modes 
during the MS (decreasing Tgff) and near the turn off (increasing 
Teff) for 13Mq models computed with overshooting. The hori- 
zontal dotted lin es stand for the observed frequency of MOST 
llSaio et al.ll2006D . The agreement between the observed frequen- 
cies and the theoretical mode spectrum is not as good as it was 
for the supergiant models. Top aov = 0.4. Bottom Qqv = 0.5. 



during the MS: overshooting, mass loss, convective criterion. 
We show that with significant mass loss or overshooting dur- 
ing the MS no ICZ appears during the post-MS phase. In 
a future work we shall extend this preliminary analysis in 
order to define an instability strip depending on those phys- 
ical aspects and compare it to the observations. In this con- 
frontation process, it is clear that asteroseismology of mas- 
sive supergiant stars can give us a better understanding on 
the physical processes not only during the supergiant phase 
but also during the MS phase. 



REFERENCES 

Aerts, C, De Cat, P., Peeters, E., et al. 1999, A&A, 343, 
872 

Alexander, D. R. & Ferguson, J. W. 1994, ApJ, 437, 879 
Bouret, J.-C, Lanz, T., & Hillier, D. J. 2005, A&A, 438, 
301 



g-modes in supergiants: mass loss and overshooting 9 



0.45 
0.40 
0.35 
0.30 
0.25 
0.20 
0.15 









15Mo 






logTe„=4.416 

\ 

\ 




■ 


\ 












. \ ^ 






v!^'" 










— 


// 

















X 






K 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
m/M 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
m/M 

Figure 11. Temperature gradients, opacity, k, and liydrogen 
abundance, X, in a model of 15Mq computed with tlie Ledoux's 
criterion for convective instability. The vertical solid lines stand 
for the boundary of the ICZ. Top MS model. Bottom Post-MS 
model. 



Martins, F., Schaerer, D., Hillier, D. J., et al. 2005, A&A, 
441, 735 

Mathias, P., Aerts, C, Briquet, M., et al. 2001, A&A, 379, 
905 

Mokiem, M. R., de Koter, A., Vink, J. S., et al. 2007, A&A, 
473, 603 

Oskinova, L. M., Hamann, W., & Feldmeier, A. 2007, 
ArXiv e-prints, 704 

Puis, J., Markova, N., & Scuderi, S. 2006, ArXiv Astro- 
physics e-prints 

Saio, H., Kuschnig, R., Gautschy, A., et al. 2006, ApJ, 650, 

nil 

Schmidt, E. G. & Carruthers, G. R. 1996, ApJS, 104, 101 
Scuflaire, R., Montalban, J., Theado, S., et al. 2007a, 
Ap&SS, 322 

Scuflaire, R., Theado, S., Montalban, J., et al. 2007b, 
Ap&SS, 431 

Unno, W., Osaki, Y., Ando, H., Saio, H., & Shibahashi, 
H. 1989, Nonradial oscillations of stars (Nonradial oscil- 
lations of stars, Tokyo: University of Tokyo Press, 1989, 
2nd ed.) 

Van Hoolst, T., Dziembowski, W. A., & Kawaler, S. D. 

1998, MNRAS, 297, 536 
Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2000, 

A&A, 362, 295 
Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, 

A&A, 369, 574 
Waelkens, C., Aerts, C., Kestens, E., Grenon, M., & Eyer, 

L. 1998, A&A, 330, 215 



Chiosi, C. & Maeder, A. 1986, ARA&A, 24, 329 
Christensen-Dalsgaard, J. & Daeppen, W. 1992, A&AR, 4, 
267 

Dupret, M.-A., De Ridder, J., De Cat, P., et al. 2003, A&A, 
398, 677 

Dziembowski, W. 1977, Acta Astronomica, 27, 95 
Dziembowski, W. A. 1971, Acta Astronomica, 21, 289 
Dziembowski, W. A., Gough, D. O., Houdek, G., & 

Sienkiewicz, R. 2001, MNRAS, 328, 601 
Iglesias, C. A. & Rogers, F. J. 1996, ApJ, 464, 943 
Klare, G. & Neckel, T. 1977, A&AS, 27, 215 
Lamers, H. J. G. L. M. & CassineUi, J. P. 1999, Intro- 
duction to Stellar Winds (Introduction to Stellar Winds, 
by Henny J. G. L. M. Lamers and Joseph P. CassineUi, 
pp. 452. ISBN 0521593980. Cambridge, UK: Cambridge 
University Press, June 1999.) 
Lebreton, Y., Montalban, J., Godart, M., et al. 2008, in 
prep. 

Lefever, K., Puis, J., & Aerts, C. 2006, Memorie deUa So- 

cieta Astronomica Italiana, 77, 135 
Lefever, K., Puis, J., & Aerts, C. 2007, A&A, 463, 1093 



